Noname manuscript No. 

(will be inserted by the editor) 



A.A. Abdujabbarov • B.J. Ahmedov • 
V.G. Kagramanova 

Particle Motion and Electromagnetic Fields 
of Rotating Compact Gravitating Objects 
with Gravitomagnetic Charge 



Received: date / Accepted: date 



Abstract The exact solution for the electromagnetic field occuring when the Kerr- 
Taub-NUT compact object is immersed (i) in an originally uniform magnetic field 
aligned along the axis of axial symmetry (ii) in dipolar magnetic field generated 
by current loop has been investigated. Effective potential of motion of charged test 
particle around Kerr-Taub-NUT gravitational source immersed in magnetic field 
with different values of external magnetic field and NUT parameter has been also 
investigated. In both cases presence of NUT parameter and magnetic field shifts 
stable circular orbits in the direction of the central gravitating object. Finally we 
find analytical solutions of Maxwell equations in the external background space- 
time of a slowly rotating magnetized NUT star. The star is considered isolated and 
in vacuum, with monopolar configuration model for the stellar magnetic field. 

Keywords General relativity • Kerr-Taub-NUT spacetime • Electromagnetic 
fields • Particle motion 



PACS 04.20.Cv; 04.40.-b; 04.70.Bw 



A. A. Abdujabbarov 

Institute of Nuclear Physics and Ulugh Beg Astronomical Institute, Astronomicheskaya 33, 
Tashkent 100052, Uzbekistan 
E-mail: abahmadjon@yahoo.com 

B. J. Ahmedov 

Institute of Nuclear Physics and Ulugh Beg Astronomical Institute Astronomicheskaya 33, 
Tashkent 100052, Uzbekistan 

The Abdus Salam International Centre for Theoretical Physics, 34014 Trieste, Italy 
E-mail: ahmedov @ astrin.uzsci.net 

V.G. Kagramanova 

Institute for Physics, Oldenburg University, 261 1 1 Oldenburg, Germany 
E-mail: kavageo@theorie.physik.uni-oldenburg.de 



2 



1 Introduction 

The existence of strong electromagnetic fields is one of the most important fea- 
tures of rotating neutron stars observed as pulsars. Starting the pioneering paper 
of Deutsch flj it was proved that due to the rotation of highly magnetized star 
the electric field is induced. The general relativistic effect of dragging of inertial 
frames is very important in pulsar magnetosphere lEUH and is considered to be a 
source of additional electric field of general relativistic origin (see, for example, 

HEME]). 

It was first shown by Ginzburg and Ozernoy [8| that an elecrically neutral 
black hole can not have an intrinsic magnetic field. However assuming that black 
hole can be inserted in the external uniform magnetic field created by nearby 
source as neutron star or magnetar, Wald [9] found for the first time the exact 
solution of the vacuum Maxwell equations for an asymptotically uniform mag- 
netic field. Then the properties of black holes immersed in external magnetic field 
were extensively studied by different authors (see, for example, |[T0lllllll2lll3|[T4l 
and for more references 1T510 . 

Despite the absence of observational evidence for the existence of gravito- 
magnetic monopole, that is of exotic space-time, called NUT space (Newman, 
Unti and Tamburino 111610 at present, it is interesting to study the electromagnetic 
fields and particle motion in NUT space with the aim to get new tool for studying 
new important general relativistic effects which are associated with nondiagonal 
components of the metric tensor and have no Newtonian analogues. 

Here we study the electromagnetic fields in the Kerr-Taub-NUT spacetime and 
in the surrounding space-time of slowly rotating magnetized relativistic NUT star. 
Our approach is based on the reasonable assumption that the metric of spacetime 
is known , i.e. neglecting the influence of the electromagnetic field on the gravita- 
tional one and finding analytical solutions of Maxwell equations on a given, fixed 
background ( estimations of contribution of electromagnetic field energy into total 
energy momentum could be found e.g. in paper 1171 ). 

We also study motion of test particle around Kerr-Taub-NUT source, which 
is immersed either (i) in an originally uniform magnetic field aligned along the 
axis of axial symmetry or (ii) in dipolar magnetic field generated by current loop. 
We use Hamilton- Jacobi equation to find influence of both NUT parameter and 
magnetic field on the effective potential of motion of test charged particles. More- 
over here we completely ignore the pathologies (existence of singularities along 
the axis or periodicity of the time coordinate to avoid them and spacetime regions 
containing closed timelike curves) of the spacetime metric due to fact that NUT 
parameter is considered as small one. 

The outline of the paper is as follows. In Section|2]we calculate the electric and 
magnetic fields generated in the Kerr-Taub-NUT spacetime following the method 
of construction of the vacuum solution to Maxwell equations in axially-symmetric 
stationary spacetime suggested by Wald |9|- 

In Sections [3] and [4] we consider the separation of variables in the Hamilton- 
Jacobi equation and derive the effective potential for the motion of charged par- 
ticles around Kerr-Taub-NUT source in a uniform and dipolar magnetic field ge- 
nerated by current loop. These results are used to obtain the basic equations go- 
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verning the region of marginal stability of the circular orbits and their associated 
energies and angular momenta. 

In Section [5] we look for stationary solutions of Maxwell equations when the 
stellar magnetic field has monopolar configuration which allows to find exact ana- 
lytical solution. Section [6] is devoted to analysis of obtained solutions for electric 
and magnetic fields exterior to the rotating compact objects with NUT parameter. 

Throughout, we use a space-like signature (—,+,+,+) and a system of units 
in which G = 1 = c (However, for those expressions with an astrophysical appli- 
cation we have written the speed of light explicitly.). Greek indices are taken to 
run from to 3 and Latin indices from 1 to 3; covariant derivatives are denoted 
with a semi-colon and partial derivatives with a comma. 



2 Kerr-Taub-NUT Source in a Uniform Magnetic Field 

We consider electromagnetic fields of compact astrophysical objects in Kerr-NUT 
spacetime which in a spherical coordinate system (ct, r, 9, ip) is described by the 
metric (see, for example, H18II191 ) 

1 2 

ds 2 = - — (A -a 2 sin 2 6)dt 2 + — [Ax - a(E + ax) sin 2 6]dtd<p 

Zj Zj 

+ ax) 2 sin 2 9 - x 2 A]dtp 2 + ^dr 2 + Sd9 2 , (1) 
where parameters £, A and x are defined by 

S = r 2 + (l+a cos 9) 2 , A = r 2 -2Mr-l 2 +a 2 , X = a sin 2 6-21 cos 9 , (2) 

I is the gravitomagnetic monopole momentum, a = J/M is the specific angular 
momentum of metric source with total mass M. 

Here we will exploit the existence in this spacetime of a timelike Killing vector 
£^ and spacelike one being responsible for stationarity and axial symmetry 
of geometry, such that they satisfy the Killing equations 

^a;fj + £,f3;u = , (3) 

and consequently the wave-like equations (in vacuum spacetime) 

□f = , (4) 

which gives a right to write the solution of vacuum Maxwell equations OA* 1 for 
the vector potential of the electromagnetic field in the Lorentz gauge in the 
simple form 

A a = dffa + CiZfc . (5) 

The constant C2 = B /2, where gravitational source is immersed in the uniform 
magnetic field B being parallel to its axis of rotation. The value of the remaining 
constant C\ can be easily calculated from the asymptotic properties of spacetime 
(HJ at the infinity. 
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Indeed in order to find the remaining constant one can use the electrical neu- 
trality of the source 



4ttQ = = i 



F af3 *dS 



a/3 



Ci 



B 



r^u a m^C {t) (uk)dS+^ f r^u a m^q v) {uk)dS 



(6) 



evaluating the value of the integral through the spherical surface at the asymptotic 
infinity. Here the equality £p. a = —£ a ;/3 = — ^Zaii following from the Killing 
equation was used, and element of an arbitrary 2-surface dS a ^ is represented in 
the form HQ 



dS al3 = ~u a A m p (uk)dS + rf^u^n^l + (ukfdS , 
and the following couples 

_ V\a^u x n>*k" VXa^u x k^m v 



k a 



y/l + {ukf ' y/1 + (uk) 2 

\uk)u a + y/l + (vJcfrf ap ''u tt mpnv 



(7) 



(8) 



are established between the triple {k, m, n} of vectors, n a is normal to 2-surface, 
space-like vector m a belongs to the given 2-surface and is orthogonal to the four- 
velocity of observer u a , a unit spacelike four-vector k a belongs to the surface 
and is orthogonal to m a , dS is invariant element of surface, A denotes the wedge 
product, » is for the dual element, ria^s is the pseudo-tensorial expression for the 
Levi-Civita symbol e Q /3 7< 5 . 

Then one can insert u = —(1 — M/r), m 1 = (1 — Al/r), and asymptotic 
values for the Christoffel symbols r? = M/r 2 and r? 3 = -3 J sin 2 0/r 2 - 1(1 - 
2M/r) cos 0/r in the flux expression (0 and get the value of constant C\ = aB. 
Parameter I does not appear in constant C\ because the integral J Q cos 6 sin 9d8 = 
vanishes. 

Finally the 4-vector potential A a of the electromagnetic field will take a form 



A 
B 



B 



X 



sm 



A 3 = -U X (a-f)+(£ + ax) 



-(S + a X )-a 2 



sm 



Using solution for vector potential $5} in spacetime fl} one could easily find 
expression for change of the electrostatic energy of the charged particle, which is 
lower to the horizon of gravitational source: 



e = eA M (^ (t) + f2\har£ M (<p))\har - eA^^ t )\ inf = elB - eaB 
At the horizon new timelike Killing vector 112 111 1411 



M + \/M 2 -a 2 + l 2 



(11) 



(12) 
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is introduced since the timelike Killing vector £ becomes spacelike inside ergore- 
gion defined by goo = 0. 

Upper limit for electric charge 



Q = 2aMB - 21MB 



(13) 



accreted by gravitational source will include in addition to the contribution from 
Faraday induction effect arising from angular momentum a (see e.g. 191. 1221 ) the 
new term being proportional to additional rotation coming from NUT parameter 
I. It is interesting to note that each parameter will accrete the charges of opposite 
sign. 

The orthonormal components of the electromagnetic fields measured by zero 
angular momentum observers (ZAMO) with the four- velocity components 



{£ + a X ) sin 2 6 - X 2 A 
AS sin 2 9 



0,0,- 



Ax - a (S + ax) sin 2 9 



AS 



(««)z 



AS sin 2 9 



V (S + axf sin 2 6 - X 2 A 
are given by expressions 

2rB sing 



(S + a X ) 2 sin 2 9 - X 2 A 
,0,0,0 



sin 9 



(14) 



E r = -- 



A 



1 - — ) S - a 2 sin 2 ' 



S^(S + a X ) 2 sm 2 9-x 2 A 
x(S + a X ) (a- |) -| [ X A~a(S + a X ) sin 2 9]^ , 



(15) 



B sin 2 9 



S 2 \jA((S + a X ) 2 sin 2 6 - X 2 A) 

{A(l + a cos 6) + 2a(S + ax - 2a 2 ) cos 9} S(S + ax) 



{S(S + ax~ 2a 2 ) - 2K(l + acos6»)} 



sin 



(16) 



B r 



Bsin6» 



SJ(S + a X ) 2 sin 2 9-x 2 A 



B" 



-{S - ax) {S + ax- 2a 2 ) cos 9 - 
2rBA 

S 2 ^A{{S + a X ) 2 sin 2 - X 2 A) 

2 ■ 2 n ( X\ 1 v-t2 -2 

I a — X ~ 2 sm ' 



xA{l + a cos 
2K 



(S + ax)(l + acos#) 



(17) 



A 



1 )S 

r 



(18) 



6 



which depend on angular momentum and NUT parameter in complex way. As- 
trophysically it is interesting to know the limiting cases of expressions ( !15l )-(fT8]l. 
for example in either linear or quadratic approximation 0(a 2 ,l 2 jr 2 ) in order to 
give physical interpretation of possible physical processes near the slowly rotating 
relativistic compact stars, where they take the following form: 

r* B Ln, *s n M(12l cos + a(l + 3 cos 29)) 1 , 1m 

E r = — i(2l + a cos 9) cos 6 —>, (19) 

r [ 2r ) 

2 Bsm6 f, / 2cos6<\ N l 
E e = \l 1 + — — +a(3cos6>- 1) L (20) 



r V sin 9 



{ 1 - ( - 4(7 2 - a 2 + a/ 



-4Z 2 cot 2 6> + a(2^ + 3a cos 0) sin (9 tan 9 j j , (21) 

B e = Bsin6»{l - — + 3__(_ a 2 + 4/2 _ 4M 2 

i r 16r 2 sin 2 9 

-8aZcos<9 + 4(7/ 2 + M 2 ) cos 2(9 + 8a^ cos36> + a 2 cos4<9) j . (22) 

In the limit of flat spacetime, i.e. for M /r — > 0, Ma/r 2 — > and Z 2 /r 2 — > 0, 
expressions d 1 5 b — (1 18b give 

lim B f = -Bcos9, lim B § = B sin , (23) 

M/r,Ma/r 2 ,l 2 /r 2 ~*0 M/r,Ma/r 2 ,l 2 /r 2 ^0 

lim £""' = lim = . (24) 

M/r,Ma/r 2 ,l 2 /r 2 ->Q M /r,Ma/r 2 ,l 2 /r 2 ^0 

As expected, expressions J23b — (124b coincide with the solutions for the homoge- 
neous magnetic field in Newtonian spacetime. 

Finally we would like to show that two form of electromagnetic field tensor 
will take more simplified form: 

B r „ T ^ „ , w „ m i n B(l + acos9) sin9 
F= — [2Kr - S{r - M)(2a - X )} w Aw + 1 

[A(S + 2a 2 - a X ) + a 2 (£ + a X - 2a 2 ) sin 2 9 - 2aK] w 2 A w° 

BrzA 1 / 2 sin (9 , „ 5 cos (9 r „ 9x 
Aw + [-«x(^ + «X - 2a 2 ) cos 0+ 

(l + a cos 0) [a(Z' + ax - 2a 2 ) sin 2 + 2zia - Z\x]] ^ 2 A (25) 



in the orthonormal Carter-type frame [22]: 

w° = (^j 1 (dt - xdifi) , w 1 = ) dr 



; 2 = S 1 ' 2 d9 , to 3 = ^jL [adt - (r 2 + a 2 + l 2 )d^\ . (26) 
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In the limiting case when NUT parameter / — » the components of field 
tensor d25b almost coincide with the expressions (3.9) of paper |9j. In the limit of 
flat spacetime, expressions for components of electromagnetic field derived using 
tensor ( 125b take Newtonian limits fl23b -d24ll. 



3 Motion of charged particles 

It is very important to investigate in detail the motion of charged particles around 
a rotating compact source with NUT parameter in an external magnetic field given 
by 4-vector potential © and d 1 Ob with the aim to find astrophysical evidence for 
the existence of gravitomagnetic charge. 
The Hamilton- Jacobi equation 

«-(£ + "'0(j£ + ■"•)-■">■ m 

for motion of the charged test particles with mass m and charge e is applicable as 
a useful computational tool only when separation of variables can be effected. 

Since Kerr-Taub-NUT spacetime admits such separation of variables (see e.g. 
QjD) we shall study the motion around source described with metric ([]]) using the 
Hamilton- Jacobi equation when the action S can be decomposed in the form 

S= -St + £tp + S ie (r,0) , (28) 

since the energy £ and the angular momentum £ of a test particle are constants of 
motion in spacetime ([]]). This is an extension of approach developed in the paper 
iflOl to the case with nonvanishing NUT parameter. 

Therefore the Hamilton- Jacobi equation ( 127b with action ( 1281 ) implies the equa- 
tion for inseparable part of the action: 

A f 9S «V ■ f d3 roV (Z + a X ) 2 sm 2 9- X 2 A ( eBK 

\dr ) + { 89 ) A S m 2 6 \ S 

2(A X -a(S + ax)sin 2 ) ( eBK \ ( e_BL\ 

zWfl { + E J V Z ) 
A -a 2 sin 2 9 / eBL\ , 

It is not possible to separate variables in this equation in general case but it 
can be done for the motion in the equatorial plane 9 = tt/2 when the equation for 
radial motion takes form 

dr \ 2 

— =£ 2 -1-2V(£, C,r,b,a,l) . (30) 
da ) 
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Here a is the proper time along the trajectory of a particle, E and C are energy 
and angular momentum per unit mass m and 



V(£, C, r, b, a, I) 



bK_\bK_ b 2 K 2 a 2 f rf_\( F bK \ 2 
+ ~MZ J M~Z + 2M 2 S 2 ~ ~£ \ + 2S ) \ ~ME J 

bK\ ( „ bL \ 



2S 



A~a 2 / bL \ 2 

+ (31) 

can be thought of as an effective potential of the radial motion which depends on 
additional dimensionless parameter 

(32) 

m 

being responsible for the relative influence of a uniform magnetic field on the 
motion of the charged particles which maybe valuable even for small values of the 
magnetic field strength ifTOl . 

Figure 1 shows radial dependence of effective potential d3~TT ) for different val- 
ues of NUT parameter I = l/M. From this dependence one can obtain modifi- 
cation of radial motion of charged particle in the equatorial plane in the presence 
of NUT parameter. As it is seen from the figure the gravitomagnetic monopole 
momentum changes the shape of effective potentials when external magnetic field 
is not strong (Fig.l, a). In the case of strong external magnetic field (Fig.l, b), 
influence of gravitomagnetic monopole momentum is negligible. 

Figure 2 shows radial dependence of effective potential ( 131b for different val- 
ues of b for fixed value of NUT parameter I = 0.5. Motion of charged particle 
in the presence of this kind of effective potential can be explained as follows: in 
the presence of external magnetic field in addition to stable circular orbits unsta- 
ble circular orbits could appear due to the appearance of maximum in the graphs 
of effective potential. From the potential we can infer the qualitative structure of 
the particles orbits. As it is seen from the figure the potential carries the repul- 
sive character. It means that the particle coming from infinity and passing by the 
source will not be captured: it will be reflected and will go to infinity again. For 
small values of electromagnetic field particles can follow both bound and unbound 
orbits depending on their energy. As external electromagnetic field increases in- 
teresting feature arises: the orbits start to be only parabolic or hyperbolic and no 
more circular or elliptical orbits exist. 



3.1 Marginally stable circular orbits 

Special interest for the accretion theory of test particles around a rotating compact 
source with NUT parameter in a magnetic field is related to study of circular orbits 
which are possible in the equatorial plane 9 = tt/2 when dr/da is zero. It is well 
known that in the presence of NUT parameter trajectory of particle lies in the cone 
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with cos 9 = ±\2l£ / C\ (see e.g. H19I ). In the case of energy and momentum of the 
particle are equal to £ = 0.9 C = 4.3 respectively, and NUT parameter I « 0.5, 
then we get 90° - 10° < 9 < 90° + 10°, which allows us to neglect deviation and 
consider the motion in the almost equatorial plane. Consequently the right hand 
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Fig. 2 Radial dependence of the effective potential of radial motion of charged particle near the 
Kerr-Taub-NUT source immersed in external uniform magnetic field for the different values of 
magnetic field in the case when NUT parameter I — 0.5 



side of equation (|3~H vanishes: 

£ 2 - 1 - 2V(£, C, r, 6, a, I) = (33) 
along with its first derivative with respect to r 

dV(S,C,r, b,a,l) 



dr 



. (34) 



In papers l23l and l24l the authors have appealed to a numerical analysis of 
analogues of equations 033b and d34K when NUT parameter is equal to zero, for 
different values of the constants of motion, orbital radii, the rotation parameter a, 
as well as the strength parameter of the magnetic field. The problem of existence 
of stable orbits in NUT space-time 11251 and other properties of particles motion 
are also discussed in our preceding paper 1261 . 

The radius of marginal stability, the associated energy and angular momen- 
tum of the circular orbits can be derived from the simultaneous solution of the 
condition 

*y( g ,Ar,M,Q =0 (35) 

or z 

and equation d33l . 

Combining equations d33l and d34l one could find the energy 

£ = -^+7±VA, (36) 
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and the angular momentum 



C = -^(l 2 + a 2 ) + i ] ±{4ab-l2 n )V\ (37) 



of a test particle, where the following notations 



2ab(r - /3) 6r - AM 

7 3r 2 -Z 2 -6/3r ' P r i2r-4M' 1 J 



A = [3r 2 - Z 2 - 6/3r] 



56 2 r 3 
4M 2 " 



(r - 4/3) 



2&V 
M 



(r - 3/3) 



3r 1- 



36 2 Z 



-2) 



(r - 2/3) 



and 



77 = 4a^7 — 6rA 



> 2 - 2/1/ 

6r7 2 - 
— 6r 



0) 



2/* 



(39) 



56 



2„3 



M 2 

36 2 Z 2 
" 4M 2 



71/ 

- 2 



2M 



(40) 



were introduced. 

Inserting now d36b and d37b into equation ( T331 ) one could obtain the basic equa- 
tion 

(r 3 - l 2 r + 2a 2 M)(j ± VA) 2 + 2a6(2/ 2 - r 2 )( 7 ± VA) 
+4aM(7 ± VA) 

x (77 ± (4a6 - 12ry)V\) + 2br 2 - l) (77 ± (4a6 - 12ry) VA) 



2^5 



6"r 



2^4 



1 



4Af 2 2M 
(a 2 - 2l 2 )r + 2M1 2 



) 



3b 2 l 2 
AM 2 J 
4aW 
M 



2M 







(41) 



The numerical solutions of this equation will determine the radii of stable cir- 
cular orbits for non-rotating NUT source immersed in uniform magnetic field as 
functions of the NUT parameter I, the angular momentum a, as well as of the influ- 
ence parameter of the magnetic field b. In the Table 1 it is shown list of numerical 
solutions for radii of stable circular orbits of particles for different values of NUT 
parameter and external magnetic field. With the increase of the gravitomagnetic 
monopole momentum, radii of stable circular orbits shifts to gravitational object, 
while external field also displace orbits to gravitational source 
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Table 1 Marginally circular orbits around non-rotating NUT source immersed in uniform mag- 
netic field 



I 







0.01 


0.1 


0.3 


0.5 


b = 


0.1 


3.66650 


3.66460 


3.66154 


3.62072 


3.53294 


b = 


0.2 


2.83292 


2.83288 


2.82876 


2.79448 


2.72020 


b = 


0.3 


2.55504 


2.55500 


2.55162 


2.52360 


2.46338 


b = 


0.4 


2.41614 


2.41610 


2.41326 


2.38984 


2,33994 


b = 


0.5 


2.33282 


2.33280 


2.33036 


2.31032 


2.26794 



4 Motion of Charged Particle in the Field of Current Carrying Loop 
Located Near the Kerr-Taub-NUT Source 



Consider a motion of test particle in the electromagnetic field created by toroidal 
currents of ionized matter rotating in accretion discs 1127 II in Kerr-Taub-NUT space- 
time. We shall employ internal solution of II28I . In particular, of the full multipole 
solution for magnetic field provided in 1128 II . we shall here focus our attention on 
the dominant (dipolar) term, which is given by 



A, 



3 fir 2 sin 2 i 



Af 3 



In 1 



2M\ 

~bTJ 



2M 



2M 



(42) 



where 



2M 



1/2 



I 



is the modulus of the dipole moment due to the current / in the loop and R is the 
radius of the loop, which is considered to be approximately equal to 6M (because 
of small values of I and a it will be not essentially modified) (see IT28II ). 

Using Hamilton- Jacobi equation d27l ) and potential d42l one can find expres- 
sion for effective potential for radial motion of charged particle in the equatorial 
plane (6 = 7r/2). It takes the following form 



eff 



1 



C V2 c, 

where we use the notations 



\ £ 2 2r l - o 



2a£r) 



c 



M 



C-2y 

2C 2 



c 



M 



(43) 



C = r 2 + Z 2 , V = Mr + l 2 , 



and 



6 en 
~8mM 2 



In 1 - 



2M 



2M 



i 



2M 



In the case, when rotational parameter a is very small effective potential takes 
following form: 



Veff=-~ 



M 



(44) 



C ' 2C 2 

which coincides with effective potential of radial motion in the equatorial plane in 
the Schwarzschild metric (see for example H29I ) if one puts £ = r 2 , r\ = Mr, 
and e = into equation (1441) . 
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Table 2 Stable orbits of test particles near the non-rotating compact object with gravitomagnetic 
monopole momentum immersed in the external electromagnetic field of current carrying loop 
(anti-Larmor orbits) 



I 







0.01 


0.1 


0.3 


0.5 


e 


= 0.1 


3.65296 


3.65305 


3.66106 


3.72475 


3.84692 


e 


= 0.2 


3.09353 


3.09360 


3.10045 


3.15506 


3.26056 


e 


= 0.3 


2.83420 


2.83427 


2.84064 


2.89154 


2.99030 


e 


= 0.5 


2.57592 


2.57598 


2.58190 


2.62931 


2.72175 


e 


= 0.7 


2.44347 


2.44353 


2.44922 


2.49488 


2.58422 



Table 3 Stable orbits of test particles near the non-rotating compact object with gravitomagnetic 
monopole momentum immersed in the external electromagnetic field of current carrying loop 
(Larmor orbits) 



I 







0.01 


0.1 


0.3 


0.5 


e 


= 0.1 


4.50043 


4.50058 


4.51456 


4.62526 


4.83569 


e 


= 0.2 


4.36234 


4.36248 


4.37677 


4.48985 


4.70439 


e 


= 0.3 


4.33040 


4.33054 


4.34495 


4.45890 


4.67491 


e 


= 0.5 


4.31295 


4.31310 


4.32758 


4.44206 


4.65896 


e 


= 0.7 


4.30800 


4.30815 


4.32265 


4.43729 


4.65445 



Plots in figure 3 show the radial dependence of effective potential, which is 
governed by equation ( 143b . for the different values of I. Like in the case of uni- 
form external magnetic field, described in the previous section, in this case grav- 
itomagnetic monopole momentum has an influence onto the effective potential, if 
magnetic field is weak (fig. 3, a). If magnetic field of loop is strong influence of 
gravitomagnetic monopole momentum is negligible (fig. 3, b). 

Consider now stable circular orbits of charged particles as it was done in Sec- 
tion[3] We shall repeat calculations, which are done in Subsection l3.ll in this case. 
Using expression (144b as effective potential we find numerical solutions for radii 
of stable circular orbits of charged particles. From results shown in the Table 2 and 
Table 3 (for anti-Larmor and Larmor orbits, respectively) one can obtain that with 
the increase of the gravitomagnetic monopole momentum radii of anti-Larmor and 
Larmor orbits shift to loop while with the increase of the electric current (creating 
dipolar magnetic field) of loop the stable orbits shift to the gravitational object. 
These results may be useful to determine NUT parameter from astrophysical ob- 
servations. 



5 External Electromagnetic Field of Slowly Rotating NUT Star for a Special 
Monopolar Magnetic Field Configuration 

In this Section we will look for stationary solutions of the Maxwell equation, i.e. 
for solutions in which we assume that the magnetic moment of the star does not 
vary in time as a result of the infinite conductivity of the stellar interior. Below 
we suggest that external electric field is generated by the magnetic field, taking as 
a special monopolar configuration. For this case we can obtain and investigate an 
analytical solution with detail consideration of the contributions from the dragging 
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1 2 3 4 5 6 

rill 

Fig. 3 Radial dependence of the effective potential of radial motion of charged particle near the 
Kerr-Taub-NUT source in the presence of current carrying loop around it for the different values 
of NUT parameter I = l/M for two cases of modulus of the dipolar magnetic field a) e = 0.3 
and b) e = 0.7. 



effects and nonvanishing NUT charge in the magnitude of the external electric 
field of the slowly rotating magnetized NUT star. 

Our main approximation is in the specific form of the background metric 
which we choose to be that of a stationary, axially symmetric system truncated 
at the first order in the angular momentum a and in gravitomagnetic monopole 
moment I. The "slow rotation metric" for exterior space-time of a rotating rela- 



15 



tivistic star with nonvanishing gravitomagnetic charge is 

ds 2 = -N 2 dt 2 + N~ 2 dr 2 + r 2 d9 2 + r 2 sin 2 ■ 

-2 [uj(r)r 2 sin 2 8 + 21N 2 cos 0] dtdtp , (45) 

that is, the Schwarzschild metric plus the Lense-Thirring and Taub-NUT terms. 
Parameter N = (1 — 2M/r) 1 ' 2 is the lapse function, u(r) = 2J/r 3 can be 
interpreted as the Lense-Thirring angular velocity of a free falling (inertial) frame. 

As a toy model we could consider the following magnetic field configura- 
tion m 

B f ~ = B f {r) ± , B § = . (46) 

Although this form of magnetic field can not be considered as a realistic, we 
will show that this toy model can be used to obtain first estimates of the influence 
of gravitational field of the NUT charge on the external electromagnetic field of 
the star. For this case, the relevant Maxwell equation reduces to 

(r*B f ) )r =0. (47) 

The solution admitted by this equation is 

B * = % > ( 48 ) 

where fj, is integration constant being responsible for the source of monopolar 
magnetic field. 

Radial magnetic field is continuous at the stellar surface and it is reasonable 
to assume that only 8 component of electric field will survive since it will be 
produced by cross product of velocity and magnetic field inside neutron star due 
to the assumed infinite conductivity of stellar matter. According to 11171 the interior 
elecric field is 

El = -e-^Bl , (49) 

where is the velocity of the stellar matter which is equal to fir sin 8 for the 
Newtonian uniformly rotating star with angular velocity S7, goo = — e 2 *. 

Then the electric field created by monopolar magnetic field is defined by the 
following Maxwell equation 



(rNE 6 ) 



„, . 2/i/cos6» f N 2 \ 
fxsinO (w ) -Jt— - — =0. (50) 
smti V ' 



The analytical solution 



a u> — Q u . 2lNcos9 a 

E e = — — ^sin0+ —-4 (51) 

civ r r sin u r z 

of equation ( 150b is responsible for the electric field of NUT star with the monopo- 
lar magnetic field d48l ). The integration constant — (Sl/i/c) sin 8 has been found 
from the matching of the exterior solution (C3 / rN) with the interior one d49b in 
the Newtonian case taking into account that the tangential components of electric 
field and the radial component of magnetic field are continuous across the surface 
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Fig. 4 Radial dependence of the ratio of the electric field to that when parameter I = for the 
different values of the NUT parameter. The value of NUT parameter I is given in cm. 

of the star. Vector potential being responsibe for the fields (I48b and d5 1 b is defined 

as 



In the figure |4] we plot the radial dependence of the ratio of the electric field to 
that when parameter I = for the different values of the NUT parameter. In this 
analysis we take typical parameters for neutron star as radius of star R = I0 6 cm, 
M = 2 X 10 5 cm, n = 2tt/(0.1s), lu = 4MR 2 f2/(5r 3 ), magnetic field at the 
stellar surface is 10 12 G. Due to the fact that in the right hand side of the expression 
d5 1 b the first and second terms have different signs the ratio is less than one. The 
results show the strong dependence of the electric field on NUT parameter. 



6 Conclusion 

Derived exact expressions d 1 5 b — <fl~8T> for electromagnetic field in the Kerr-Taub- 
NUT spacetime indicate that electromagnetic field will be affected by the gravit- 
omagnetic charge. However the induced electric field dl5l ). dl6l ) depends on NUT 
parameter I linearly while the magnetic field ( 1171 ), d 1 8b depends on I quadratically. 

Analytic general relativistic expressions for the electromagnetic fields external 
to a slowly-rotating magnetized neutron star with nonvanishing gravitomagnetic 
charge I are presented. The star is considered isolated and in vacuum, and for 
simplicity with the monopolar magnetic field directed along the radial coordinate. 

We have shown that the general relativistic corrections due to the dragging of 
reference frames and gravitomagnetic charge are not present in the form of the 
magnetic fields similar to dipolar case I8l l31l but emerge only in the form of the 
electric fields. In particular, we have shown that the frame-dragging and gravito- 
magnetic charge provide an additional induced electric field which is analogous 
to the one introduced by the rotation of the star in the flat spacetime limit |6). 




(52) 
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Motion of charged particles around Kerr-Taub-NUT source immersed in ei- 
ther (i) external uniform or (ii) dipolar magnetic field have been investigated using 
Hamilton- Jacobi equation. We have shown that in the presence of NUT parame- 
ter and magnetic field the shape of effective potential will be changed. However 
modifications caused by external electromagnetic field are dominating. Investiga- 
tion of the stability of motion of charged particles shows, that external magnetic 
field shifts orbits of test particles to gravitational source in both cases, while NUT 
parameter shifts to gravitational source in the case of uniform magnetic field and 
towards loop (in opposite direction) in the case of the presence of current loop 
around the gravitational source. 
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